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(Dated: ) 

In recent years Bayesian Inference has become an efficient tool to analyse single molecule trajec¬ 
tories. Recently, high density single molecule tagging, Langevin Equation modelling and Bayesian 
Inference [lQ] have been used to infer diffusion, force and potential fields at the full cell scale. In 
this short comment, we point out pitfalls Q,i to avoid in single molecule analysis in order to get 
unbiased results and reliable fields at various scales. 

PACS numbers: 


A response to Holcman et al jjj] 

Understanding the physical and biological processes governing the trafficking of membrane proteins is an cen¬ 
tral question in cell biology and biophysics. A powerful approach to model the motion of a random walker in 
a complex environment is to use an overdamped Langevin equation. In the case of a coarse-grained description 
(corresponding to the limited spatial and temporal resolution accessible in experiments), this equation writes: 

r = F (r) /q (r) + V 2D ( r h (*) 

where F (r) /y (r) corresponds to a local drift, 7 (r) is a spatially-varying friction, D (r) a spatially-varying diffusion 
coefficient and rj (t) is a zero-averaged gaussian noise. 


In past years, the advent of single molecule techniques has opened new possibilities to estimate the values of the 
local parameters F(r),7(r) and D (r) in the cellular context. In particular, bayesian inference tools have been 
developed to extract local force and diffusivity fields based on single molecule trajectories iH. Not only these tools 
have been extensively characterized from a methodological point of view Sill but they have also been applied to 
experimental systems, in order to investigate the dynamics of toxin receptors in lipid rafts 0 , @, fj] and to calibrate 
optical tweezers (7]. 


More recently, the analysis of the force and diffusivity fields has been extended to the full cellular scale in [ 2 ] 
and in thanks to high-density single molecule techniques, such as sptPALM [12] and uPAINT flji]. Thereby, 
it has been possible to quantitatively analyze the synaptic stabilization of AMPA receptors ( 2 ] and Glycine 
receptors M- Both papers similarly interpreted the role of synaptic scaffolds in terms of local confining potentials in 
line with a view of the neuronal membrane where neuroreceptors diffuse and transiently stabilize at synaptic sites .Hl- 


In their comment, Holcman et al. question the validity of the approach and description presented in 0 - We reply 
below to the points that they raised by: (i) addressing some computational issues (statistical estimation of the model 
parameters, numerical simulations) and demonstrating biases in the statistical estimators used in [ 2 J, (ii) discussing 
the interpretation of the results. 


Parameter estimation and numerical simulations. 

A key point for the determination of the force and diffusivity fields based on single molecule data is the choice of 
the statistical estimators. To clearly illustrate the simplifying assumptions used in [2j, we first remind the standard 
definition of the simplest likelihood used to analyze the dynamics of individual biomolecules in the absence of positional 
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corresponding to N translocations Ari (during the sampling time At) in a subdomain characterised by a constant 
drift f = F /7 and diffusion coefficient D. By differentiating log(L) with respect to f and D, one directly gets: 


and 


<91og(L) 

df 


N 


E 


(Ari — fAt) 
2D 


dlog(L) 

dD 


N 

~D 


N 

E 


(Ari ^ fAt ) 2 
AD 2 At 


The maxima are obtained for: 


N 

f = (^ An)/NAt 

i=l 


and 

* j 

The estimators used in Q correspond to the maxima in the limit fAt <C Ar, and can thus be designated as 
Maximum Likelihood Estimators (MLE). 

The estimators used in Q are biased when applied to experimentally measured trajectories for two main reasons. 
First, an important element in the analysis of single molecule trajectories is the positioning noise cr, which results 
from a combination of processes such as photon noise, motion blur, electronic noise and limited performance of peak 
fitting algorithm. However, the effect of positioning noise was entirely neglected in the analysis presented in 2], 
which, as shown below, leads to estimator biases and to improper parameter values for the synapses studied in 2 ]. 
Second, the approximation fAt <C Ari (not discussed in |2|) also leads to strong bias in the force and diffusion 
estimation over range of values relevant for the experiments reported in j2|. In our work, we have used Bayesian 
inference schemes and directly included in our computation the full positioning noise as an experimental parameter 
a (see eq. equations 5 and 6 in 0)- We refer to Q for a demonstration of the optimality of Bayesian infer¬ 
ence methods for diffusion measurements and to [5;] for a detailed discussion on their use in the case of confined motions. 

To directly compare the performances of the estimators used in Q and 0 , we have performed numerical simulations 
corresponding to experimental conditions reported in [1] (a = 25nm, At = 20ms, about 40 points per subdomain). 
We first estimated the diffusion coefficient in the absence of external drift for D varying between 0.01 and 0.5 /im 2 /s 
(Figure 1A). For D < 0.05/im 2 /s (a value comparable to the diffusion coefficient in synaptic domains), we observe a 
clear bias when using the MLE estimators whereas our estimator remains accurate over the entire range. The limit 
corresponds to the regime where a is comparable or greater than V DAt. 

Next, we tested the estimation of the diffusion coefficient and drift when brownian particles (with a diffusion 
coefficient D = 0.05/zm 2 /s) were submitted to a force varying between 10 ~ 2 and 10 1 pN. While our approach yielded 
accurate estimates of F and D , there was a systematic bias in the MLE estimation of both the force and the diffusion 
coefficient at low force. Furthermore, a strong deviation can be observed for forces larger than F c ~ 0.1 pN (Figure 
l.Bl,B2), a value corresponding to the the limit where the displacement due to the drift fAt. is comparable to the 
diffusional movement V DAt. Importantly, F c is comparable to the confinement force Uq/25 in a potential well with 
depth Uo = 8 kgT and extension <5 = 200 nm, such as the ones reported in Overall, our simulations indicate that 
the estimates of the parameters (diffusion, forces and potential) for the synaptic domains are not accurate in [ 2 } 

Note that Bayesian inference techniques have been recently applied to evaluate the experimental trapping force 
in optical tweezers [7| - This system is advantageous since the trap parameters can be finely tuned with the light 
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intensity and, importantly, independently calibrated. The estimate of the trap spring constant determined with 
inference schemes was in good agreement with (in fact more accurate than) the results obtained by standard tools 
(equipartition and power-spectrum methods). 

Another important point is the extraction of the potential from the force field. In Q|, this is done by minimizing 
the square difference between the gradient of the potential and the force field. Such class of optimization problems, 
especially with a discrete set, limited number of points and harmonic approximation, are always solved with a 
regularization term [161. Il7l| (see equation 7 in [lOj). Solving these problems without any regularization, either of 
the force field or of the potential, leads to biases (as illustrated in Figure 1C,D). 

Finally, regarding the numerical simulations (including the Gillespie scheme), we refer to [l4j for approximating the 
Fokker-Planck by the master equation and to [15f | for simulating the master equation. We especially emphasize the 
importance of the Gillespie scheme for rapidly generating a large number of trajectories, which enabled us to compute 
fundamental descriptors of the motion (such as the propagator) as function of the time and over a large range of 
distance. 

Interpretation and discussion of the model. 

A legitimate topic of discussion is the interpretation of the diffusivity and force (or potential) fields estimated 
from individual trajectories. In their experiment^ Hoze et al noted the existence of domains where individual 
translocations pointed toward a single center point [2[. Similarly to what has been reported earlier for toxin receptors 
ii, such behavior was described in terms of a trapping potential. For AMPA receptors, the potential wells were 
interpreted as resulting from interactions of the receptors with scaffolding proteins PSD-95 via stargazin molecules. 

Before discussing and comparing our results 0 , we emphasize that our inference approach does not require 
the computation of a diffusion and a potential ( D,U ) but can be performed to estimate the parameters (D,F) 
or (D^F/'y) as it was done in The inference scheme simply provides a description of the environment in 

terms of the diffusivity and force fields that have the highest probability. In the specific case of our experiments 
on glycine receptors 10], we have favored a description in terms of a potential for several reasons. First, the force 
field did not exhibit rotational terms, thus justifying the computation of a potential. Second, glycine receptors are 
known to interact with gephyrin scaffolding proteins at inhibitory synapses. The fact that the trapping areas of the 
membrane proteins perfectly co-localized with the location of gephyrin clusters strongly supported the notion that 
the clusters acted as trapping wells. The role of the scaffolds in the observed potentials was further demonstrated 
experiments using membrane proteins with mutant constructs of the /3-loop mediating the receptor-gephyrin 
interactions: for a mutant with weaker affinity, we measured shallower traps, and for proteins lacking the /3-loop 
entirely, there was a complete absence of potentials. In this context, the potential depth represents the energy needed 
for receptors to escape gephyrin clusters. It is worth noting that it is the classical definition of the potential in physics. 

In floj , we have used a spatially varying friction 7 (r) , reflecting the heterogeneity of the plasma membrane at the 
full cell scale. Furthermore, we imposed the relationship D (r) = ksT/'y (r). This hypothesis was supported by the 
fact that the transitions between adjacent domains (determined experimentally) satisfied detailed balanced condition, 
thus suggesting a local equilibrium. For more advanced discussions on the FD relationship for systems with spatially 
varying friction, we refer to recent references 

An important difference between the results of Q and floj | lies in the quantitative measurements of the potential 
depths. For the wild-type beta-loop, we found an average energy of 3.4 k#T with less than 15% of the traps above 6 
k^T (measured over ~ 69 clusters). In contrast, 25 % of the traps reported in [3] were deeper than 8 k^T, indicating 
a much higher level of stabilization. While it could be due to the difference between the biological systems and their 
molecular constituents, it seems likely that the biased estimators (see discussion above) and the limited number of 
data points presented in @ (2, 5 and 10 wells depending on the experimental conditions) also contribute to the 
discrepancy. Additional experiments, combined with careful statistical analysis, are needed to clarify the difference 
between excitatory and inhibitory scaffolds. 

In our view, the main remaining challenge is to quantitatively determine how the trapping energies at scaffolding 
clusters relate to the biochemical properties of their individual constituents. We emphasize that there is no general 
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FIG. 1: Comparison of the results of the estimators used in 2| and 0 in simulated trajectories. All values are averaged over 
1000 trajectories. The theoretical values are in black, the results associated to the estimators described in [gj are in red, and to 
our estimators in cyan. Positioning noise was set to 25 nm f2j. A Evolution of the estimated diffusion with constant diffusion 
field matching experimental conditions in [2|. Note that blue and red error bars are the same size, distortion is due to log-log 
plot. B1,B2 evolution of the force and diffusivity (respectively) with force for constant diffusion field D = 0.05;im 2 .s _1 . C 
Example of an extracted potential from a simple ID force field. The mesh was 50nm 1- D evolution of the difference between 
the input potential and the reconstructed one for 2D force fields, as a function of the anharmonicity. 


procedure to go from a diffusion and force fields to biochemical characteristics of complex molecular assemblies. It will 
depend on many geometrical and molecular properties of the assemblies, as well as on the spatio-temporal resolution 
of the measurements. Fortunately, nanoscopv imaging methods now provide powerful tools to analyze the composition 
and structure of scaffolding complexes 0 El and relate them to their functional roles at synapses. 


Conclusion. 

In summary, we stand by the methods and the approach described in 0 and earlier papers, and point out to errors 
in the statistical and computation methods used in [I}, [2| . Note also that we did not comment here on other important 
elements for the mapping of the diffusivity and force (or potential fields), such as the meshing procedure and the 
mesh size which should be adjusted to the averaged translocation length (a point not discussed in 0), the effect of 
strong variation between the number of points between neighboring mesh domains, the effect of confinement on the 
accuracy of parameters extraction, the effect of heterogeneous diffusion field on the convergence of the estimators 
and the projection of the potential on a proper functional basis. We refer the reader to previous publications 
(including their supplementary materials) for further discussion on all these points. 
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